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We propose a model Hamiltonian for charge transfer along the DNA double helix with temperature 
driven fluctuations in the base pair positions acting as the rate limiting factor for charge transfer 
between neighboring base pairs. We compare the predictions of the model with the recent work of 
J.K. Barton and A.H. Zewail (Proc. Natl. Acad. Sci. USA, 96, 6014 (1999)) on the unusual two-stage 
charge transfer of DNA. 



PACS numbers: 87.15.-v, 73.50.-h,82.30.Fi 

Charge transport and electrical conduction are known 
to occur in a wide range of organic linear chain crys- 
tals of stacked planar molecules Jlj]. Transfer rates 
from molecule to molecule are determined by the single- 
particle integral r, with typical rates of the order of 10 15 
sec -1 . Strong interaction between the electronic degrees 
of freedom and molecular vibrations may reduce this to 
10 12 sec -1 , a typical lattice or intramolecular vibration 
frequency. By comparison, biochemical charge transfer 
processes, such as those encountered in the metabolic re- 
dox (oxydation- reduction) chains 0, usually are much 
slower (down to 1 sec). Key steps often involve some 
form of large scale motion of the molecule. 

DNA can be considered as a one-dimensional, a- 
periodic, linear chain of stacked base-pairs. More than 
30 years ago it was suggested Q that duplex DNA might 
support electron transport in a manner similar to that of 
linear chain compounds, namely by tunneling along over- 
lapping 7r orbitals located on the base-pairs. Barton et 
al. Q first presented evidence that photo-induced, rad- 
ical cations can travel along DNA molecules in aqueous 
solution over quite considerable distances (more than 40 
A, or about ten base-pairs). If so, DNA might present 
us with flexible, molecular-size wires able to transport 
charge in aqueous environments. Possible applications 
range from micro-electronics to long range detection of 
DNA damage. Subsequent studies by a number of groups 
reported a wide range of values for the effective inverse 
spatial carrier decay length f3, ranging from as low |^] 
as 0.02 A -1 to as high as around one A -1 , the differ- 
ent values most probably reflecting differences in charge 
transfer for different base sequences 0. 

Recently, Barton and Zewail [|| (BZ) used femtosecond 
spectroscopy to measure the rates of the DNA charge 
transfer process. An unusual two-step decay process was 
observed with characteristic time scales of 5 and 75 ps 
respectively. Ab-initio molecular-orbital calculations || 
find that DNA has a large single particle band-gap, and 
a transfer integral r of order 0.1 eV. This would lead to 
a charge transfer rate r/h for coherent tunneling that is 
comparable to that of the linear chain compounds, but 
that is much too high compared to the rates measured 
by BZ. 

Apart from coherent tunneling B, a number of al- 



ternative transport mechanisms have been proposed, in 
particular incoherent, phonon-assisted electron hopping 
between bases, with the electron wave fully localized on 
each subsequent base-pair ■ This would reduce the 

transfer rate to a typical intra-molecular vibrational fre- 
quency (ps -1 ), but this is still much too large to explain 
the slow second stage step of the decay. It was also sug- 
gested H that a charged radical could induce a polaronic 
distortion Jl^] of the DNA internal structure that might 
control charge transfer. The explanation proposed by BZ 
for the long relaxation times is that large amplitude ther- 
mal fluctuations of the intercalated photoreceptor sets 
the rate limiting step for the charge transfer. 

The aim of this letter is to construct a model Hamil- 
tonian to treat charge transfer along a chain under con- 
ditions of large structural fluctuations, and suggest that 
thermally induced structural disorder interferes with the 
it orbital overlap mediated charge transfer, leading to 
long relaxation times. To construct this Hamiltonian, we 
first must discuss the origin of the structural fluctuations. 
Figure 1 shows an example of a typical DNA configura- 
tion obtained by a molecular dynamics (MD) simulation 
]i~3[ . The relative orientation of neighboring bases along 
DNA is characterized by a set of collective variables such 
as the relative roll and twist angles (R and T), and the 
relative slide displacement (see Fig.l). Long time MD 
simulations of DNA lead to typical RMS fluctuation an- 
gles for R and T of order 5 and 9 degrees jl4} in the ps to 
ns time-window, while the mean base-pair spacing also 
shows large amplitude fluctuations Structural fluc- 
tuations in the ps to ns time-window have been observed 
|l6f experimentally as dynamic Stokes shifts in the fluo- 
rescence spectrum of the DNA. The local fluctuations are 
extraordinarily strong compared with those due to ther- 
mally excited phonon modes in crystalline linear chain 
materials. The unusual 'softness' of the R,S, and T vari- 
ables is also reflected by the fact that their mean values 
vary greatly depending on base-pair sequence Q . 

In our simplified model, we include only two collective 
modes. The first mode is an angular variable 9(t), which 
is that relative rotation angle of the two bases which 
couples most efficiently to the 7r-orbital tunnel matrix 
element. Next, the displacement variable y(t) represents 
that collective mode which couples most efficiently to the 
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FIG. 1. Top: Snapshot of the result of a 140 ps room 
temperature MD simulation of DNA, by Swaminathan (left 
panel) as compared to the ordered DNA structure at t = 
(right panel). Bottom: Schematization of the roll and twist 
angles R and T and of the relative slide displacement variable 
S. 

on-site energy of the radical. This second form of cou- 
pling provides the necessary mechanism for energy trans- 
fer between the charge and the the thermal reservoir re- 
quired for hopping transport along a random sequence 
of base-pairs with different on-site energies. Both 0(t) 
and y(t) are treated as classical harmonic variables that 
are coupled to a heat bath of oscillators. As a result, 9(t) 
and y(t) obey -in the absence of the radical- the following 
Langevin equations: 

My + ly y + Mtf y y = r) y (t), 

i9 + le e + m 2 g 9 = m (t). (i) 

In Eq. ^, / is the reduced moment of inertia for the rela- 
tive rotation of the two adjacent bases, £lg is the oscillator 
frequency of the rotation mode, while M and £l y are the 
reduced mass and natural frequency of the displacement 
mode. The values of M,I, Qg,Q y , and the damping co- 
efficients are obtained by comparing the Fourier power 
spectra of y{t) and 6(t) obtained from Eq. |l|, to power 
spectra of MD simulations of DNA pi . From a typical 
long time (10 ns) MD series, we find oscillation periods 
of 27r/f2 of order 1-10 ps, a (large) mass M of order 1- 
10 kDalton, (1 Dalton = 1 a.m.u.) a moment of inertia 
/ of order 10 2 ks^lg 2 , and relaxation times comparable 
to the oscillation period (i.e. the slow modes are close 



to critical damping). The amplitudes of the white-noise 
variables r\ (t) follow from the fluctuation-dissipation the- 
orem for classical variables. We will neglect mode cou- 
pling between different pairs of adjacent base-pairs. The 
two modes are then coupled to a one-dimensional, tight- 
binding Hamiltonian for single particle charge transport: 

H = ^{r(6» i! i + i)(c l t +1 c i + c\c i+ i) + ejcjct 

i 

+ ~M i (y J 2 + ( w + yoA^) 2 )} + H bath ({9, y}). (2) 



In Eq. g, €i is the on-site electronic energy. The distance 
?/o,i is the change in the equilibrium value of the y variable 
of the i-th base when the particle localizes on that site, 
while Mft 2 yQ is the typical deformation energy. 

Certain limiting cases of this general Hamiltonian are 
familiar from studies of one-dimensional charge trans- 
port. For uniform e.j and for fixed 9, H is the Hamil- 
tonian of a tight-binding polaron ||l2|| . For fixed 9 and 
and y and random , H is the Anderson Hamiltonian for 
localization in one-dimension. For the case of DNA, we 
assume that site-to-site differences in the value of are 
of order O.leV based on the sequence dependent differ- 
ences in the ionization potential |19|. Next, the transfer 
integral t(6) will be assumed to be small compared to the 
thermal energy k B T for 9 near a special value, denoted 
by 9*, the 'rapid decay state'. Finally, the characteristic 
interaction energy MVl 2 y^ between the charged radical 
and the on-site structural variable y is assumed to be 
large compared to fcgT (e.g. due to electrostatic effects) 
and of order e. 

In this unusual, high temperature, strong coupling 
regime, the transfer integral t(9) is the lowest energy 
scale. Under these conditions, particle motion described 
by H is indeed dominated by incoherent hopping from 
site to site. We first restrict ourselves to the case of a par- 
ticle which resides at site A at time t — and then hops 
to the neighboring site B with a different on-site energy. 
For fixed 9, the transition rate T(9) for incoherent charge 
transfer between A and B can be computed by applying 
the method of Garg, Onuchic and Ambegaokar (GOA) 
pof . In Fig. 2 we show the two potential energy surfaces 
V+(y) = (A\H\A) and V~(y) = (B\H\B) for the particle 
respectively on the A and B sites as a function of y. Effi- 
cient transfer between the two potential energy surfaces 
takes place nearly exclusively at the crossing points y* 
where V + (y*) — V~(y*), shown in Fig. 2. Note that an 
energy barrier E r must be overcome to reach this cross- 
ing point. In the high temperature, strong coupling limit, 
the on-site probability decays exponentially with a rate: 
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FIG. 2. Potential energy surfaces for a particle localized 
on either of two neighboring sites A and B as a function of 
the collective variable y. The 9 variable is assumed fixed. 
Non-adiabatic transfer between the two surfaces takes place 
exclusively at the degeneracy point y* . 



where the energy scale E r defined in Fig. 2, depends on 
the difference between the on-site energies. The validity 
condition for Eq. || is that t(9) must be small compared 
to (EVfcsT) 1 / 2 . By itself, Eq. || does not account for a 
two stage decay, as observed by BZ. In a typical ensem- 
ble of radical sites (with the same on-site energies) there 
must be significant heterogeneity, with less likely states 
characterized by low charge transfer rates. This hetero- 
geneity is incorporated by demanding that t(9) is appre- 
ciable only for 9 near the special value 9* , while tunneling 
plays no role for different 9 values 



21 1 . If 9 undergoes 



large amplitude thermal fluctuations, then there should 
be considerable heterogeneity for the transfer rates. 

Assume then that, at time t = 0, an ensemble of par- 
ticles is prepared on the A site, with the 9(0) variable 
obeying the Boltzmann distribution. We define the prob- 
ability density P(9, t)d9 to be the fraction of radicals at 
time t that are still on site A, and whose 9 values is in 
the range between 9 and 9 + (19. For the overdamped case 
of Eq. 0, P(9, t) obeys the following equation: 



dP{0,t) k B Td 2 P(9,t) 1 d 



at 
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d9 2 
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[m 2 e 9P(9,t)} 

-T(6)P(6,t). (4) 



We can now use Eq. [| to discuss the decay rate. For 
F(9*) large compared to the thermal equilibration rate 
Tg 1 two forms of decay are encountered p3j. 

(i) Early stage decay: At times t = 0, a certain frac- 
tion of the oscillators has an energy exceeding E — 
(1/2)1 (fl g 9*) 2 . These oscillators will pass through the 
9 = 9* point within a time of order max{rg, Qg 1 }- When 



this happens, there is a finite probability for charge trans- 
fer to take place. On a time scale of order max{re, Qg 1 }, 
these high energy oscillators are removed from the prob- 
ability distribuition. 

(ii) Late stage decay: After the high energy oscillators 
have been removed from the distribution, further decay 
requires energy diffusion along the oscillator scale from 
lower energies towards E(9*) = (l/2)I(il s 9*) 2 . Once the 
energy of an oscillator reaches this value, efficient charge 
transfer takes place. After a standard, but lenghty, anal- 
ysis of Eq. [| we find that the late stage decay rate is: 



kLATE — Tgilg 



(5) 



The factor in front of Eq. [s] is of the same order as the 
early stage decay rate while the term in brackets is of the 
order of the thermal probability that 9 exceeds 9*. The 
second stage decay rate strongly increases with increas- 
ing temperature, while the early stage decay rate is not 
expected to be strongly temperature dependent although 
the fraction of sites that exhibits early stage decay should 
be strongly temperature dependent (of the same order of 
the term in brackets of Eq. |5|) . 

The time dependence of P(9, t) would be consistent 
with the observations of BZ if this thermal probability of 
the 9* state is of order 10~ 2 . In that case, one percent 
of the sites would show rapid decay with time scales of 
order ps, while the remainder would show decay slowed 
down by a factor of 10 2 . We only treated here the nearest 
neighbor hopping process. Charge transport over longer 
distances described by our Hamiltonian reduces - under 
the assumed conditions - to a classical one-dimensional 
diffusion in a random medium with site specific transfer 
rates. The transport properties of such systems have 
been extensively discussed elsewhere |P3fl . 

In conclusion, we propose that charge transport along 
DNA proceeds by classical diffusion with high amplitude 
thermal fluctuations providing the rate limiting step for 
the site-to-site charge transfer. If correct, charge trans- 
port along DNA would have unique characteristics as 
compared to linear chain compounds. Since the radical 
severly deforms the local structure it might be considered 
as a polaron in the strong coupling limit but, unlike pola- 
ronic transport, hopping is controlled by thermal fluctua- 
tions. Indeed, Eq. || predicts that the charge transfer rate 
should strongly increase with temperature which is consi- 
tent with the observations of BZ. A better description of 
the mechanism proposed in this paper for charge trans- 
fer along DNA would be to consider it as a repeated se- 
quence of reversible oxydation-reduction reactions. The 
site-to-site charge transfer would be viewed as a chemical 
reaction dominated by a transition state where the col- 
lective variables y and 9 assume a special value (y* and 
9* respectively). We are not aware of any of the linear 
chain compounds exhibiting this curious form of charge 
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transfer. On the other hand, a recent single-molecule op- 
tical study of a particular reversible oxydation-reduction 
reaction did report pi[ two-stage non-exponential behav- 
ior but with decay rates much lower than those measured 
by BZ (in the range of 1 sec -1 ). 

The higher rates of charge transfer in DNA would be 
due to the fact that the molecular motion of the bases still 
is significantly restrained by the backbone. If the present 
analysis is appropriate, then charge transport in DNA 
occupies a unique position intermediate between charge 
transport in solid-state materials and charge transport in 
biochemical charge transfer reactions. 

Finally, the proposed Hamiltonian obviously incorpo- 
rates a number of rather serious simplifications. We in- 
clude the collective modes only in a schematic way. We 
require a large structural on-site distortion of a site by 
the particle (of order 0.1 eV), but this is likely to require 
an anharmonic description of the collective modes. We 
did not include coupling of modes of adjacent pairs, the 
double-stranded nature of DNA with the possibility of 
inter-chain charge transfer or effects related to the ter- 
tiary structure, i.e. the coiling of the duplex. 
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